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Stationary, axisymmetric and asymptotically flat spacetimes of dust of which trajectories are 
integral curves of the time translation Killing vector are investigated. The flow has no Newtonian 
limit. Asymptotic flatness implies the existence of singularities of the curvature scalar that are 
distributions and that are not isolated from regularity regions of the flow. The singularities are 
closely related to the presence of additional stresses that contribute negative active mass to the 
total (Komar) mass, which is zero for asymptotically flat spacetimes. Several families of solutions 
are constructed. 

o : 

O , I. INTRODUCTION 

<N ■ 

We study stationary, axisymmetric and asymptotically flat spacetimes of dust in free fall along integral curves of 
D ' the time translation Killing vector. For reasons which will become clear later, we shall call this motion 'van Stockum 
Mh , flow'. Despite the flow is a rigid rotation with zero angular velocity with respect to asymptotically static observers, 
00 ■ the squared vorticity scalar (which for the particular flow equals, up to a constant factor, the proper energy density 
CNJ \ of dust) in contrast to Newtonian physics, does not vanish. Angular velocity of the flow with respect to locally 
non-rotating observers numerically equals the angular velocity of dragging of inertial frames, while angular velocity of 
matter in linearized gravity is many orders of magnitude greater. This shows that van Stockum flow is ultra-relativistic 
even in the limit of negligible density. 

We demonstrate that asymptotic flatness implies the existence of curvature singularities that have distributional 
character, spatial measure zero, and are not isolated from regularity regions. The singularities are closely related 
£C) • to additional weird stresses with negative active mass. Total mass of such spacetimes is necessarily zero, which has 
already been conjectured by Bonnor 0] , and total angular momentum is determined by the amplitude of the dipole 
in a multipolc expansion. 

We construct two classes of asymptotically flat solutions and a class of asymptotically non-flat solutions. Bonnor's 
C"| . solutions belong to the first two, and the van Stockum solution V\ to the third. 

Or 

6 ■ 
a: 

We shall focus on axisymmetric, stationary and asymptotically flat spacetimes of dust flowing along opened integral 
curves of the time translation Killing vector By asymptotic flatness the axial symmetry Killing vector rj, of 
which integral lines are closed, vanishes on the symmetry axis at least for radii sufficiently large. Consequently, i) 
V[a£.(3^n;^] = and £[ a r)pr) tt . 1 ,] = at least at one point. The energy-momentum tensor is proportional to hence 

ii) ^Rf}" £, a ri^ = and r] fJ 'R t }' J £ i a r]^ = on the basis of Einstein's equations. In addition, we assume Hi) [£, rj\ = 0. 
Under these three assumptions a theorem proved in Q guarantees that, maybe apart from isolated points (p, z), there 
exist a coordinate system in which the line element takes the general form 

ds 2 = -V (di - Kdc/)) 2 + V" Vd^ 2 + e 2 * (dp 2 + Adz 2 ) , (2.1) 

where V, K, and A are structure functions of two variables p and z. In these coordinates £ and r\ attain the 
particularly simple form £ M = <5f and rf = 5%. 

The four-velocity of van Stockum flow reads u — Z£, thus Z~ 2 = —^S,^ and £ M Z p = 0. Killing equations imply 
u (ii;v) = £(pZ,v)i hence the expansion scalar u^.^ = 0. On projecting u^-v) onto the it-orthogonal subspace and taking 
the traceless part, one infers the shear tensor vanishes identically, as well, van Stockum flow is therefore rigid. Despite 
the fact and that angular velocity of the flow also vanishes (tx oc £), the vorticity scalar u> 2 does not, where 

2 _ 1 (^) 2 (VS) 2 q ^ 
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II. VAN STOCKUM FLOW 



Typeset by REVTeX 



2 



A. Equations of van Stockum flow 

By definition of the flow, the energy- momentum tensor reads 7j,„ = DZ 2 ^^, where D is the proper energy 
density such that £ M 2? „ — = rfT) ^. Einstein's equations and the contracted Bianchi identity imply the flow 
is continuous and geodesic. As = 0, the continuity equation (VZ^) — is satisfied identically. In addition, 
u v u KV = —u u Z£ U]IJ/ . The geodesic equation u v u KV = will be satisfied if = —u v u v . p + u u u u (YnZ) . As u^u p = —I, 
Z must be constant. This, in turn, implies V is also constant and, without loss of generality, we may set V = 1. On 
defining = g-rfrf + + rf^£ v we obtain for dust (T M „ - Tg^/2) = 0, then Einstein's equations 

imply K^R^ = or pe~ 2 *<9 p In y/\A\ — in coordinates, hence A = A(z). If so, the form of (|2.1[) allows us to set 
A(z) = 1. We have thus shown that the line element of van Stockum flow reads 

ds 2 = -cl< 2 + 2K (p, z)dtd<f> + (p 2 - K 2 (p, z))A<p 2 + e 2 *< p ^' (dp 2 + dz 2 ) . (2.2) 

Let E ii v = W v - \R^ V - 8?rT^, then E p p = and E p z = yield the following relations 

K 2 , - K 2 K n K , 

•* - -"V- •- - —rf- (2 - 3) 

The integrability condition \P „ z — ^ , zp imposes on K the elliptic constraint 

CK = 0, C = d 2 --d p + d 2 . (2.4) 
P 

Provided (|2 .3(1 and (|2.4|) are satisfied, the other components of E p, v , but E t t and E*^, vanish identically. The latter 
two will also vanish for C 2 solutions (we stress, the reservation 'for C 2 solutions' is necessary) if only 



K 2 + K 2 
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III. ASYMPTOTIC FLATNESS AND CURVATURE SINGULARITIES 

A function K e is a regularized profile of a solution K of (|2.4|) in an open subset V of the plane (p, z) if K e £ C°°(V) 
and K e — > K as e — > almost everywhere in V. We shall denote by X the set of all irregularity points of van Stockum 
flow, that is, a subset of M 3 where CK does not exist for solutions in the sense that lim e ^o f x dmfCK e ^ for 
any S > 0, where X C X s , < dist(dX 5 ,X) < 6, e 2 */ = p^K c or p 2 K e and dm = e 2 ^ pdpdtfidz . In particular, 
K <^ C 2 (X). We recall the result of the theory of elliptic equations that X/S 1 has measure zero in the plane (p, z). 

In what follows we shall prove that 1^0 for asymptotically flat van-Stockum spacetimes. For if we suppose for 
contradiction that X = 0, then inside a ball Br C M 3 bounded by a two-sphere Sr of radius R and concentric with 
the origin 

De 2 *pdpAd0Adz £ = = ° — / — (K z dp — K „dz) A d<^> = — / — r— 7-d0 A d6 (3.1) 

on J p 8n J smv 

Br Sr Sr 

in virtue of the Stokes theorem, provided (K 2 )^ r — o(sin#), (rsvciO = p, rcosr? = z). By asymptotic flatness 
K ~ 2 Jr _1 sin 2 ^asr-> oo, hence, for R sufficiently large, the surface integral on the right-hand side is negative and 
tends to as R — > oo, while the volume integral is positive, a contradiction, thus indeed X ^ 0. 

The surface integral (|3.1[) coincides in the limit R — > oo with the total (Komar) mass M, which is zero. Indeed, for 
asymptotically flat spacetimes with metric (|2.2j) . M is given by 1 

M = -— lim / -y^ £Q V ^rd^ Q A dx r3 = lim — / -^--^ A d0 = °- ( 3 - 2 ) 

Sr Sr 



1 The 2-forms in 113.21 1 are equal on <Sr irrespective of whether Einstein's equations are satisfied or not 
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We conclude, therefore, that the integral liniR^oo J g dm(8np 2 ) x e 2 ^KCK, which was omitted deliberately in (|3.ip 
as it would be zero for C 2 solutions, does not vanish. For regularized profiles K e this integral tends in the limit e — > 
to minus the total mass j R3 ^ x T>dm of the regularity region ]R 3 \ X. Putting this in other words, asymptotically flat 
van Stockum spacetimes contain additional sources of negative active mass located in X. The sources are distributions 
as they have both finite mass and a spatial support of measure zero. 

To be more explicit, the proper energy density T> — T^u^u" , the trace of spatial stresses S — T^„ (u ,1 u v + g^"), 
Tolman's active mass density on a hypersurface of constant time Dt — {8tt)~ 1 and the curvature scalar 

R = R^^ of a spacetime with metric (|2.2p read 

V = | e -2*£l|±^ _ i e -2» + * „) , V T = e -2*Ek±£k + i e -™*CK, 



For any smooth ^ and X satisfying only (|2.3p . the above definitions reduce to 

P = P + e- 2 * I ^/:^, V T = V + e-™^CK, 

lbTrp 7 p 

In particular, for regularized profiles K e we obtain, outside Is and in the limit e — ► 0, I? = 2?t = D = (87r) 1 i? 
and iS = 0, like for dust. However, in the same limit J T Sdm ^ for arbitrarily small 6. Note, that X is the scalar 
curvature singularity as R is a distribution on X. Indeed, if K e is a regularized profile of an asymptotically flat solution 
then J R3 R ^ / H 3\j R = 8ir j R3 V in the limit e — > 0. Thus i? = 87rP + jx where D is smooth and integrable, and 71 

is a distribution localized on X. The singularity is not isolated from regularity regions. 2 

Here we give an example illustrating the statements. The Bonnor solution [3j with an embedded surface 
layer of negative mass, and which we shall denote by Kg, can be regularized by defining K^(p,z) — ^/8a 3 p ■ 

p 2 ■ ((a + Vz 2 + e 2 ) 2 + p 2 y 3/2 , a > 0, fi > 0. Although K% is globally C°°, its limit K% is not even dif- 
ferentiable in X, which is the plane z = 0. On integrating over R 3 and taking the limit e — > we obtain 
J R3 Ddm = f R3 ^ x Ddm — /i, J R3 Ddm = 3/i/4, Jj Sdm = /i/4, J R3 Rdm = Airp and J R3 Dxdm = 0. The lat- 
ter holds identically as %-KyJ—gDTdp A d<fr A dz = d [jjT x KK , z dp A d<fi + p^ 1 KK p d(p A dz) for C 2 functions, hence 
J R3 DTdm = M = for regularized asymptotically flat profiles. Since for Ifjj 8ttD = R only outside T, and since 
J^ij-Ddm = J R3 Ddm = /1 ^ /i/2 = (8tt) 1 J R3 Rdm, the curvature scalar is a distribution, and R is smooth and 
bounded only outside X. 

These properties and those discussed in the introduction show that global and asymptotically flat van Stockum 
spacetimes are not viable physically. The same, for example, concerns also asymptotically flat stationary and ax- 
isymmetric spacetimes of which internal metrics are matched onto the external line element (|2.2p . as then M = 0. 
However, a possibility that van Stockum spacetime can be part of a regular spacetime cannot be excluded. 

IV. THREE CLASSES OF SOLUTIONS 

A. Solutions with a layer of negative mass 

Solutions to equation (|2.4p can be sought via integral transforms, for example, 

00 00 
K{p,z) = p [ \Kj(X)e' xlz ^ 1 (\p)d\, or p [ XK K (X) cos (Xz)K 1 (Xp)dX 



generate z-symmetric solutions; Ji and Ki are Bessel functions. To give an example, the solution p 2 r 3 discussed in 
[2j] has Kj(X) — 1 and Kk(X) = 2/7T, while the Bonnor solution K% has Kj(X) oc e~ aA . The latter belongs to a class 



2 In this sense we cannot agree with the statement of paper Q that a solution found therein had no curvature singularity. 



of solutions defined by specifying Kj{\) — l ^ n +i)\ e aA > 71 <= which yields 

K „ (p , 2)= ( ? ±i^V. 2Fl f3 +nj2 + „ ;2i ^ 



(a+ |z|) 2 »+ 3 * 1 \2 ' ' ' (a+\z\Y 

Apart from the plane z = the solutions are smooth everywhere. One can show that \K n \ < p(l/ (a + |z|)) 2 ™ +2 , 
thus, at least for I < a, hypersurfaces of constant t arc globally space-like as then \K n \ < p. Since K n d r K n (r,8) ~ 
r -(4n+3) s j n 4 q ^i mes a bounded geometrical factor, the spacetimes are asymptotically flat with M = 0. Function 
De 2 * is finite for z ^ and for r sufficiently large behaves as r ^( 4ll + 6 ) . The plane z — is thus a curvature singularity 
with finite and negative active mass. Only the Bonnor solution Kg (n = 0) has nonzero angular momentum. 



B. External and Internal multipolar solutions and a multipole expansion 



Let K(p,z) — W (r)Y (cos 9) , where p — rsin#, 



rcosO. There exist three families of solutions to equation 



satisfying r 2 W"(r) = \W(r) and (1 - x 2 )Y"{x) + \Y(x) = 0, (x = 



cost 



— cos 2 (a)/4, < a < it/2 and Hi) 
internal (W = r 



with i) A = a(a + 1), a > 0; ii) 



cosh (a)/4, a > 0. The i) class contains x-analytic external (W 



1 ) solutions. In this way we obtain external Kg and internal Kj multipolar solutions 



l ) and 



0,2,4,6, 



1,3,5,7,. 



(p 2 +z 2 Y 



(p 2 +z 2 ) 



—1 



z(p 2 + z 2 y /2 A„ 



U + z>) {n+1)/2 B n 



where A n {y) = 2 F X (± + f, -§; §, y) and B n (y) = 2 F X 



— n • — z 



Internal solutions Kj^ (of which 



2 ' 2 ' p 2 +z 2 

element is the van Stockum solution give rise to spacetimes that are not asymptotically flat. With the exception 
of the monopole Kg , K^ yield asymptotically flat spacetimes that contain in the center pathological singularities 
with non-integrable V (e.g. e 2 *2? > a 4 / (2-Kr 6 ) for the dipole Kg = a 2 sin 2 (0)r _1 ). For the solutions M = 0, thus 

contributions to M from the singularities are formally — oo. Another non-physical property of K^ is that \K (p, z)\ > p 
in the vicinity of the center. In the region the axial symmetry Killing vector rj is time-like, as so, the region contains 
closed time- like curves (e.g. Kg is such inside the region bounded by p — a\/sma sin a, z = a cos aVsina, a S (0, it)). 

However, external multipoles Kg" 1 ^^ appear in multipolar expansions of asymptotically flat solutions. We illustrate 
this by giving an example below. 

It should be clear that K a — — J" a a~ 1 sdsf(p, z — s), a > 0, where /(p, z) = KgP = zr~ l , is a z-symmetric solution 
such that < K a < a. The conformal mapping z + ip = a cosh (u + iv) (u > 0, < v < n) is invertible apart from 
two points (p, z) = (0, ±a) and transforms S a = {(p, z) : p = 0, z G [—a, a]} to the segment u = 0, v G [0, tt]. In this 
map the solution reads 



if Q (u,v) = asin 2 (u) ( cosh(u) + - sinh 2 (u) In 



tanh z ( - 
2 



The resulting spacetime is asymptotically flat as K a ~ (4/3)ae "sin u for large u. In the vicinity of u = 
_ftT a ~ u 2 lnu, therefore X = S a . Multipolar expansion of K a in the base of functions K^\ n = 1, 3, 5, ... , reads 

2 p 2 2 1 p 2 (p 2 - 4z 2 ) 4 3 p 2 (p 4 - 12p 2 z 2 + 8z 4 ) 6 



K a(p, z) = -—a 
6 r J 
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and * a (r, 6») - -(a 4 /72)r -4 (7 + 9cos26>) sin 2 6 as r -> oo. Asymptotically, ([2~2]) reduces to 

4„2 

ds 2 - -dt 2 sin 2 OAtAcj) + dr 2 + r 2 (d6 2 + sin 2 (9d</> 2 ) . 

3 r 

Comparison with the asymptotic expansion of the Kerr metric gives total mass M = and total angular momentum 
J = a 2 /3 in agreement with (|3.2[) and with the analogous expression for the total angular momentum of asymptotically 
flat van Stockum flow 
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